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Abstract. - We study the non-equilibrium dynamics of the spherical p-spin models in the 
scaling regime near the plateau and derive the corresponding scaling functions for the corre- 
lators. Our main result is that the matching between different time regimes fixes the aging 
function in the aging regime to h(t) = exp(i 1_M ). The exponent fi is related to the one giving 
the length of the plateau. Interestingly 1 — /i is quickly very small when one goes away from the 
dynamic transition temperature in the glassy phase. This gives new light on the interpretation 
of experiments and simulations where simple aging was found to be a reasonable but not perfect 
approximation, which could be attributed to the existence of a small but non-zero stretching 
exponent. 



Introduction. - Glassy phases are widely observed in Nature and in materials [1]. The 
most stricking features are very slow relaxations and history dependent phenomena, such as 
rejuvenation or memory [2,3]. Such phenomena reflect that the system is aging: the older 
the system, the slower the relaxation. The arrest of the dynamics is generally attributed to 
large scale freezing, and the approach to the glass transition has recently been shown to be 
associated with strong dynamic heterogeneity [4]. Some regions of space may relax very fast, 
while others may relax very slowly. In order to relax almost frozen regions, many particles 
have to move cooperatively, and the dynamical arrest arises when the time for relaxing large 
regions - the relaxation time - exceeds the available experimental time. The understanding 
of the freezing phenomenon and of the low temperature - or glassy - phase has been boosted 
by the introduction and the resolution of models such as the spherical p-spin model, or the 
mode-coupling theory [5]. In these models, the dynamics stops at some critical temperature 
T c , the idealised glass transition temperature. The dynamics of the spherical p-spin model 
has been solved in the limit of very large waiting times [6], which lead to the following 
phenomenology in the glassy phase. In this phase a sharp separation of timescales occurs, with 
two distinct time sectors [7] , which in turn correspond to fast and slow degrees of freedom. The 
former correspond to local (in time) equilibration of local (in space) degrees of freedom and 
is often pictured by a rattling inside "cages" , which respond to a small external perturbating 
field through Fluctuation Dissipation Theorem (FDT). The latter rule the slow structural 
relaxation, following a slow descent in the energy landscape [8], and their response can be 
described using a modified FDT, in which the temperature T is replaced by an effective 
temperature T* > T . Still, old samples exhibit stationarity on finite time differences. Another 
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striking consequence of the analytic solution of the p-spin model in the limit of an infinite 
waiting time is that for very old samples the correlation function C(t, t') (we will always choose 

t > t' , with t written on the left) takes the form C (jj^tfj when t — » oo with t ~ if . This 

form has been used to fit data in experiments and simulations of aging systems. In many 
cases, h(t) is fairly well approximated by a power law ("simple aging") [9, 10] ; however better 
fits have often been obtained by using h(t) = exp (£ 1-M /(1 — fi)), with /j, < 1 but fj, close 
to 1 ("sub aging") [9,11,12]. The interpretation was that the times used for measurements 
were two small and that longer waiting times would provide better agreement with simple 
aging. However, the study of the p-spin (or related) models has not helped resolving this 
issue, as in this case h(t) is left unknown by the dynamic equations in the aging regime. In 
other words, in the aging regime, the dynamics is invariant under strictly increasing time 
reparametrizations [6,13]. 

In this Letter, we revisit the analytic solution of [6], by unveiling a new time scale, which 
as we shall show is crucial in order to understand the crossover from the stationary to the 
aging regime, and which we shall call the "scaling regime" . It corresponds to times of the 
order where t w is the waiting time and f3 a new exponent, and was studied in detail in the 
spherical Sherrington-Kirkpatrick (SK) model in [17]. Our main results are the calculation of 
the aging function h(t) and the description of the large but finite waiting time behaviour of 
the dynamic correlators. 

The model. - We focus on the celebrated spherical p-spin model is defined by the Hamil- 
tonian H p = ^2 1<il ... i <N Ji 1 ---i p 4 > h ' ' ' fapi where TV is the number of soft spins, which verify 
the spherical constraint J^i 4*1 — N and the J^.-.j 's are Gaussian random coupling with van- 
ishing average and variance N p /p\. Except when specified, we shall consider only p > 2, 
where a discontinuous ideal glass transition occurs [14]. The dynamic equations for these 
models in the large N limit involve the correlator C(t,t') = (<f)i(t)<j>i(t')) and the response 
R(t,t') = ( s%'(t') ) ^° an ex ternal field rji- It is convenient to use the integrated response 
F(t, t') = L ds R{t, s), which leads to the following equations [7]: 

(d t +z(t))F(t,t') = l + J dsE(t,s)F(s,t'), (1) 

(d t + z(t))C(t,t') = f dsY,(t,s)C(t',s)+ f dsD(t,s)d s F(t',s), (2) 
Jo Jo 

where E(i, s) = -&£=&C(t, s)P- 2 d s F(t, s), D(t, s) = §C(i, s) p " 1 and the Lagrange multiplier 

z(t) = T + ^ J^ > dsC{t,s)P~ 1 d s F{t,s) enforces the constraint J2 i M t ) 2 = N - Th ere is no 
characteristic timescale in Eqs. (|H2p . as the time is expressed in units of the microscopic 
timescale t$. 

The different time regimes. - In the limit of large waiting time t! — + oo, the dynamics can 
be partly solved assuming the sharp separation between the stationary and aging regimes [7] . 
In this limit, when t—t' ~ 1, the solution of (JTJ) and © is stationary: C(t, t') — Csrit — t') and 
R(t,f) = R ST (t - t'), where C ST and R ST are related by FDT: TR st (t) = -dC ST {T)/dT 
(t > 0). The stationary part of the correlation behaves at large r as Cst{t) = qEA + 

c< ST T ~ a + c< ST T ~ 2a + - "> wh ere a verifies ^-2a) = T/(2T*). On the other hand, when 

t - t' ~ if, the solution has the form [6] C(t,t') = C (j$f), F(t,t') = -^(^y), which 

satisfy !F(X) = The effective temperature T* verifies T/T* = (p ~ 2)( ^~ 9ea) , while 
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the Edwards-Anderson parameter satisfies 2 Q%a(^ — Qea) 2 = T 2 . Near A = 1: C(A) = 

qEA — c^ G (l — X) b — c^(l — A) 26 H , where b = 1 for the spherical p(> 2)-spin model (we 

keep b when possible, as some of the results apply to a wider class of models) [15]. 

Now we add an extra timescale to the picture. We assume the existence of an exponent f3 
which governs the size of the plateau with t' , or equivalently t. We describe time differences of 
order t@ , interpolating between the stationary and aging regimes, thus (3 < 1. The existence 
of this scaling regime is confirmed by numerics [16] and by the analytic solution in the case 
p = 2 [17]. The thickness of the plateau is described by another exponent a such that C(t, t') = 
q E A+t- a giC-=f) + t- 2a g 2 ( t -ff) + ---aiidF(t,t') = ^A+t^ a w 1 { t -^f) + t- 2a w 2 ( t -^f) + --: 

This implies R(t,f) = t^r^-^f) + i~ 7 ~ Q r 2 (^) H , with 7 = a + f3, and r 4 (x) = 

dwi(x)/dx. In addition, one also introduces T(x) — —g[(x)/w' 1 (x), which goes from T at 
small x to T* at large x ; the finiteness of T(x) explains the occurrence of the same a for C 
and F. 

Matching the different regimes. - The scaling regime interpolates between the stationary 
regime at small arguments and the aging regime at large arguments. This gives relations 
between different exponents and scaling laws. First we take r = t — t' = et 13 , with t — > 00 and 
e <C 1 (in this order). Matching with the stationary regime gives t~~ a gi(e) ~ e _>o c"g T e~ la t l P a , 
and thus <?i(e) ~ Cg T €~ ia , a — [3a. It also gives io,(e) ~ —Cg T /Te~ la . This is similar to 
the analysis around the plateau in the high temperature phase [5], Second, we take r = xt@ 
with x 3> 1 and match with the aging regime. On one hand C(t,t') — C ( J — °ba — 

c^(T</)(t)) h , with 0(t) = dln ^ (t) . On the other hand C{t,t') = q EA + t- a gi (x) + ■■■. 

Comparing the dependences in t and t gives gi{x) ~ — c^g.(>lia:)'' i , Wi(x) ~ c^Q(^4ia;) h Yr* 
and <p(t) ~ t _AI /Aj, with /i = /3 + a/6. On can conclude that /i(f) ~ exp ^xTCl" 7 /!) ) ' P rov ided 
4>(t) — t~^ <C 1/t. The corrections to are given by the corrections to the leading order 
of the Wi's (i > 2) [18] and can be neglected when /i + a > 1. There is numerical evidence 
(see below) that this is the case for p-spin models, and thus cf>(t) = t~^/Ai + o(i _1 ). This 

gives h(t) ~ exp ^xf\r^77))' anc ^ thus < fj, < 1, as /i(t) must map infinity onto infinity. 

We set Ai to 1, as changing A\ to A\ is equivalent to changing C(A) to C(A j4l Z' 41 ). This 
freedom to choose A\ reflects that the miscroscopic timescale is forgotten in the aging regime. 
It provides further support for the absence of non-vanishing (except logarithmic) corrections 

to \nh{t) - ^-^/(l-fi). Indeed, ]nh(t) = i 1_A1 /(l - /i) + A 2 t 5 H , with S > would fix the 

microscopic timescale to through A 2 oc tg . This argument can even be used to give directly 
the form of h(t). Indeed, assuming C(A) = A^ to) , h(t) = exp/(t/t ) and d to C(t, t') = in the 
aging regime gives f(u) — A(t n ) + B(t )u c [18]. Thus h is either a stretched exponential or 
a power law, while other empirical forms [12] are ruled out. Remark that in order to obtain 
h(t) by matching the different regimes, we have integrated the large time behaviour, so this 
in principle gives h(t) only approximately. However the meaning of the aging function is only 
asymptotic, and valid for large t. Note that the relations we have derived between exponents 
are very general, provided the expansions we used for the late stationary and early aging 
regimes are valid. 

Equations for the scaling regime. Now we sketch the - lengthy - derivation of the 
equations for the scaling functions at the plateau, which will be detailed elsewhere [18]. The 
following procedure gives equations for gi and u>i, but the exponents a and /3 go away, and 
P remains unknown. We set t — t' = xt 13 in |T|) and ^ and separate the integrals in different 
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time sectors. Then, the identification of different orders gives the equations. In the mean time, 
one gets for large times z(t) = 2(00) — z\t~ a — Z2t~ 2a — • ••, where z\ = and z% > [18]. 
The multiplier z(t) being a one-time quantity, its value is rather easy to get in numerics, 
and thus provides a useful estimate for a ; in addition, it also gives the behaviour of the 

energy density £{t) = (T — z{t))/p. Now we write / Q * ds = J* + f*~ M + + X t '_ e ^' 

which separates t — s in several distinct regimes: very short waiting times, aging, scaling and 
stationary. Next we take the limit t — » 00, followed by rj — > 0, A — > 00 and e — » 0. We 
then identify all orders in t~ a : taking the above limits leads to divergences which cancel each 
other giving en passant the values of the coefficients of the leading orders of the g^s and w^s 

(i > 2). We also write J t , — J t , +Clt + $ t ,+ * t p + f t _ e t p and send t\ and e 2 to zero. We now 
give the results of the identification of different orders. First, the orders t° and t~ a of (|1I2|) 
give: z(oo) — p ^~ 1 ^ g|T^ 2 (l — qEA) + 1 q ^ A , and we get at order t~ 2a : 

w(x) 2 + f dyw'(y)^+ f dyw'(y)(w(x-y)-w(x))=C F , (3) 

Jx 1 JO 

dy [(p - l)g{y) {w'(x +y)- w'(y)) + w'(y) (g(x +y)- g{y))] 
+ (p-2) I dyw'(x + y)(g(x + y)-g(y))+ f dyw'(y){g(x-y) + {p-2)g{y)) (4) 

Jx 

dy [(p - l)g(y)w'(x + y) + w'(y) (g(x + y) + g(x -y)- 2g(x))} 

+ + 2w(a:o)5W y* — V EA 9™ X + 2T* = C ' 

where xq is such that g(xo) — 0, and g^, Cf and Cc will be given somewhere else [18] (y is 
given below). In order to obtain consistent equations, we had to make the assumption that 
C{t,t') vanishes faster than t~~ 2a when t' is finite and t — > 00. This assumption is consistent 
with the stretched exponential behaviour of hit), and has been checked numerically. We thus 
ended up with coupled equations for 171 and w\. It can be easily checked that (J3j is compatible 
with the asymptotic behaviours obtained from the matching with the other regimes ; the 
same check for Q is harder. In addition, integrating these equations from small to large x 
in principle provides the value of The p-spin model is a special case where the form of 
the aging scaling function is known ; indeed one has simply [6]: C(A) = qEA^ u , which gives 
c^q = —v (recall A = 1). Thus, C is in principle determined by g\ and w\, although we have 
not yet been able to solve (|3|4|) analytically or numerically. 

The case p — 2. - For p = 2, C(t,t') and R(t,t') can be computed exactly [19], and 
T* = 00. In addition the exponents of the scaling regime are [17] j3 = 4/5 and a = 2/5. 
Moreover the analysis must be slightly refined [18], as g\ = and correspondly c}q = ; it 
gives fj, = 1, which corresponds to simple aging, as confirmed by the exact solution. Remark 
that here one has Z2 = 0, while Z3 — —3/8. 

Numerical check. As one exponent is still unknown, we have checked the relations 
between different exponents using a code for numerical integration of the equations provided 
by the authors of [20]. We start with p — 3. At T = 0.5, we obtained a « 0.3 from the fit 
of z(t), and (3 ~ 0.66 from the scaling of tp{t') which solves C(t' + tp{t'),t') — qEA- This is 
in very good agreement with a = (3a (a rj 0.448) and gives /1 RJ 0.96, which provides a very 
good collapse of the data in the aging regime, as found empirically in [16]. The same test was 
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Fig. 1 - Exponents a (circles), f3 (triangles) and /i (squares) versus temperature for p — 3 (solid lines), 
p — 4 (dashed lines) and p — 10 (dotted lines) . 



carried out at several temperatures from T = to T = 0.61 (T c = y/3/8 ~ 0.6125), providing 

very good agreement with our analytic description. In Fig. [TJ the exponents a, (3 and \x are 
shown (along with the same for p = 4 and p = 10). The observed tendency is that they 
decrease with p and T. In addition, j3 and \i saturate respectively to 2/3 and 1, their values 
at T = 0. The value /z(T = 0) = 1 corresponds to approximate aging slightly above T = 0. It 
is clear from Fig. [T] that /i + a > 1 is a reliable working assumption. Finally, we have checked 
numerically the collapse of the data when using the scaling we introduced for the plateau 
regime: in FigO —TdF/dC is represented versus (C(t' + t, t') — qEA)t' a parameterised by t 
for several waiting times, for p = 3 and T — 0.5. A consequence of our scaling hypothesis is 
that a master curve is obtained, identical to the curve of T/T{x) versus g(x), which shows 
the onset of FDT violation during the scaling regime. 

Discussion. - In many glassy systems, fits of the aging regime were done using h(t) ~ t 
(or any power law), giving reasonable agreement. It was however observed that h(t) ~ 
exp gives better collapse of the thermo-remanent magnetisation in spin-glass exper- 

iments, with fi slightly below 1 [9,11]. An interpretation of this fact was that the waiting time 
was not large enough for the aging scaling to be seen well. In the case of the p-spin model, 
as we have shown, the stretched exponential form is related to the power law dependance of 
the timescale of the plateau with the waiting time and to the forgetting of the microscopic 
timescale in the aging regime. In addition, we have found \i < 1 from numerical estimates, 
which appears as an evidence against simple aging. However /i departs significantly from 1 in 
the vicinity of the critical point, and the numerical estimation of the exponents may be flawed 
by critical effects. For instance, the exponents shown on Fig. [T] seem to have a cusp at T C) 
which could be an indication that the scaling regime is strongly influenced by the cusp of qEA 
at T c [5], as the position of the plateau itself fluctuates a lot at the approach to T c . A simple 
way to discriminate such effects would be to compute the variation of (3 with T. Indeed our 
numerical results suggest that p decreases with T, while (3{T) obtained by imposing fi = 1 
at all temperatures - which gives (3{T) = 1/(1 + a(T)) - is an increasing function. It remains 
that for fixed t', our data and those of [16] are clearly compatible with stretched exponen- 
tials rather than power laws. The scenario fj, < 1 is at least relevant for the time regimes 
numerically accessible. 
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Fig. 2 - Test of the scaling hypothesis for p — 3 and T = 0.5. Inset: same as main, but with 
C(t' + T,t') as horizontal coordinate ; the crossing occurs at q and corresponds to r = xot 13 , with 
g(x ) = 0. Waiting times are t' = 2197, t' = 4295, t' = 8590, t' = 17180 and t' = 34359. 

Conclusion. - In this Letter, we have shown that the matching of the stationary and 
aging regimes in glassy systems requires the precise description of the intermediate scaling 
regime, which is defined at large but finite waiting time. We have shown that the natural 
assumption that in this regime the correlators obey power law scaling leads to a stretched 
exponential form for the aging function. In addition, the exponent is related to the one giving 
the duration of the plateau. The assumptions and ansatze which have been made in order to 
get information about the scaling regime have been checked with very good accuracy using 
numerical integration of the exact dynamic equations. The consequences of these results are 
threefold. First, it shows the p-spin model as an archetype system where the aging function 
is close to a power law, but different, at least in the accessible time regimes. This provides a 
motivation for revisiting experimental and numerical data fits where the same phenomenon 
was observed. Second, our calculation provides explicit corrections at large but finite times to 
the solution given in [6] . This may indicate the route to follow in order to study the dynamics 
of systems where the ideal transition does not really exist, and where an aging behaviour, 
with a waiting time dependent effective temperature was found. Here a calculation at large, 
but not infinite waiting time is clearly necessary. Third, we hope that our treatment may lead 
to a proper mathematical approach to the solution of the dynamic equations. Indeed, it has 
recently been proven rigorously [21] that the dynamic equations governing the dynamics in 
the (thermodynamic) limit N — > oo derived in [6] . The important challenge is now to compute 
the exponent (3, i.e. the length of the plateau. Indeed, it is important to confirm whether or 
not n < 1 is not a numerical effect due to critical fluctuations. Finally, our analysis remains to 
be adapted to models with continuous replica symmetry breaking such as the SK model [22]. 
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